Energy thresholds for discrete breathers 
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Discrete breatliers are time-periodic, spatially localized solutions of the equations of motion for 
a system of classical degrees of freedom interacting on a lattice. An important issue, not only from 
a theoretical point of view but also for their experimental detection, are their energy properties. 
We considerably enlarge the scenario of possible energy properties presented by Flach, Kladko, 
and MacKay [Phys. Rev. Lett. 78, 1207 (1997)]. Breather energies have a positive lower bound if 
the lattice dimension is greater than or equal to a certain critical value dc. We show that dc can 
generically be greater than two for a large class of Hamiltonian systems. Furthermore, examples 
are provided for systems where discrete breathers exist but do not emerge from the bifurcation of 
a band edge plane wave. Some of these systems support breathers of arbitrarily low energy in any 
spatial dimension. 

PACS numbers: 45.05.+X, 63.20.Pw, 05.45.-a 



The phenomenon of locahzation is of interest in almost 
any branch of physics. In addition to the well known An- 
derson localization due to disorder, the last decade has 
seen an increasing attention for localization phenomena 
in translational invariant systems, i.e., in the absence of 
disorder. Discrete breathers are time-periodic, spatially 
localized solutions of the equations of motion for a sys- 
tem of classical degrees of freedom interacting on a lat- 
tice. A necessary condition for their existence is the non- 
linearity of the equations of motion of the system, and 
the existence of discrete breathers has been proved rig- 
orously for some classes of systems 1', '3, Q . In contrast 
to their analogues in continuous systems, the existence 
of discrete breathers is a generic phenomenon, which ac- 
counts for considerable interest in these objects from a 
physical point of view. In fact, recent experiments could 
demonstrate the existence of discrete breathers in vari- 
ous systems such as coupled optical waveguides 'i^ , low- 
dimensional crystals antiferromagnetic materials 6], 
Josephson junction arrays (0 and references therein), 
and molecular crystals [g]- 

Properties of discrete breathers, as well as of some gen- 
eralizations of discrete breathers, have been studied in de- 
tail in a large variety of different models. However, apart 
from existence proofs and studies of the spatial localiza- 
tion of discrete breathers (which is typically exponential), 
only a few general results exist. Among these it is worth 
mentioning the remarkable result by Flach, Kladko, and 
MacKay |9j on energy thresholds for discrete breathers 
in one-, two-, and three-dimensional lattices (and subse- 
quent generalizations to systems with long range inter- 
actions (l^ and to partially isochronous potentials pj]|). 
These results on energy thresholds have practical rele- 
vance, as they can assist in choosing a proper energy 
range for the detection of discrete breathers in real exper- 
iments or computer experiments. The main achievement 
of these papers on energy thresholds can be summarized 
as follows: for a certain class of Hamiltonian systems of 



infinite size, a critical spatial dimension dc exists, such 
that 

• for a system whose spatial dimension d is smaller 
than dc, no positive lower bound on the energy of 
discrete breathers exist, i.e., discrete breathers of 
arbitrarily low energy can be found, 

• for a system where d > dc, there exists a positive 
lower bound on the energy of discrete breathers. 

The critical dimension dc depends on the type of non- 
linearity present in the system and is, according to 
typically two, but never greater than two. This result 
was obtained by tacitly assuming that the Hamiltonian 
admits a Taylor series expansion around its equilibrium 
point. Although this is a reasonable assumption for many 
systems or models in physics, there exist some remark- 
able exceptions (see e.g. and 01). In the present 
Letter, some assumptions made in 9] are critically re- 
visited and a larger class of Hamiltonian systems, not 
necessarily expandable around the equilibrium point, is 
considered. We find a much richer scenario for the exis- 
tence of energy thresholds for discrete breathers: 

1. The critical dimension dc can be greater than two 
for large classes of Hamiltonian systems. We report 
on systems which do not show an energy threshold 
for discrete breathers in spatial dimensions two and 
three. 

2. There exist discrete breathers which do not emerge 
from a tangent bifurcation of the band edge plane 
wave. For these discrete breathers an energy 
threshold exists which, however, cannot be ob- 
tained by the analysis used in and 0, 0, Q| . 

3. In the absence of a linear spectrum ("sonic vac- 
uum"), discrete breathers are superexponentially 
localized and dr — oo. 
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We consider a d-dimensional hypercubic lattice with 
N sites. Each site is labeled by a c?-dimensional vector 
i G li'^, and a state Xi E M.^ is assigned to each lattice 
site, where / is the number of components and is to be 
finite. The Hamiltonian of the system is of the general 
form 

= [^loc(^0 + ^int(a;i, {xi+j})] , (1) 

where Tij^t describes the interaction of a state at site i 
with the states {2;^+^} in a (finite) neighborhood. Ti. is 
assumed to have an equilibrium point at Xi — 0, with 
^({0}) 0- III order to exemplify, in what follows we 
shall make reference to the form 



p1 



V{q^ 



, (2) 



where is a local potential, W an interaction potential, 
and Ni denotes the set of nearest neighboring sites of 
i. V and W are assumed to be two times continuously 
differentiable (C^) and are allowed to have an isolated 
non-analyticity at zero, being of the shape 



V{q) = V2q' 



n— 3 



(3) 



where the Vn are constants which may be different for 
q>0 {V+) and g < (V^). (Analogous for W.) 

For generic Hamiltonian systems, periodic orbits — and 
hence also discrete breathers — occur in one-parameter 
families, and some typical choices of quantities to index 
such a family are the energy E, the amplitude measured 
at the site with maximum amplitude, or the breather 
frequency cjf,. It may or may not be the case that, in a 
certain system, discrete breathers with arbitrarily small 
amplitude can be found. (This is in contrast to 0, where 
it is argued that the breather amplitude can always be 
lowered to arbitrarily small values.) 

Systems with small amplitude breathers and phonons: 
First, it is important to recall that, for discrete breathers 
to exist as generic solutions, the discrete breather fre- 
quency LUb has to be nonresonant with the phonon spec- 
trum ujq, i.e., ojq/iOb ^ N Therefore, in the limit of 
small amplitudes, the discrete breather frequency u)b bas 
to approach an edge of thephonon spectrum from out- 
side the phonon spectrum j^]. With amplitudes tending 
to zero, it is observed that localized solutions (discrete 
breathers) and extended ones (phonons) become more 
and more similar, and one might suspect that discrete 
breathers appear through a bifurcation from a phonon 
mode. For some one-dimensional systems this has been 
proved 0], in other cases numerics supports this conjec- 
ture. Calculating the critical energy Ec at which this bi- 
furcation occurs, the minimal energy of discrete breathers 



in a system can be obtained. Flach computed the 
critical energy for systems where the potentials V and 
W are infinitely many times differentiable {C°°) around 
their minima. Extending this type of analysis to sys- 
tems where the Hamiltonian is in the vicinity of its 
equilibrium point, one obtains 



Er cx N^-^/'^ 



(4) 



relating the bifurcation energy Ec to the number of lat- 
tice sites N . Note that, as C°° is a subspace of C^, also 
the result by Flach, Kladko, and MacKay [Ec oc iV^-^/d-j 
can occur, however, as a non-generic case (not robust un- 
der C^-perturbations of the Hamiltonian). Additionally 
to the above mentioned physical applications of nonan- 
alytic potentials, the fact that the generic behavior is 
somewhat ruled by these potentials renders our result 
interesting from a theoretical point of view. 

Considering the proportionality Q in the limit of an 
infinite number of lattice sites, 



lim Ec 



=0 if <i<4, 
>0 if d>4. 



(5) 



different energy properties are obtained, depending on 
the spatial dimension d of the system. For d < 4, dis- 
crete breathers of arbitrarily small energies can be found. 
For c? > 4, however, breather energies do not approach 
zero even if the oscillation amplitude goes to zero, and 
below a certain positive energy threshold no breathers 
exist. The value of the threshold energy can be obtained 
from perturbation theory following the ideas of So, 
extending the result of Flach, Kladko, and MacKay j^, 
for this larger class of Hamiltonian systems we obtain 
a critical dimension dc = 4, instead of dc = 2, below 
which generically there is no positive lower bound on the 
discrete breather energy. 

The mechanism which leads to the existence or non- 
existence of an energy threshold can be understood some- 
what intuitively in the following way: when lowering the 
amplitude of a discrete breather towards zero, two differ- 
ent mechanisms, competing with respect to their effect on 
the breather energy, take place: the localization strength 
of the discrete breather tends to become weaker, leading 
to an increase in energy, while the mere lowering in am- 
plitude causes a decrease in energy [see Fig. ^a) for an 
example] . Depending on the respective strengths of these 
effects (which in turn depend on the spatial dimensional- 
ity of the system), an energy threshold may exist or not. 

To confirm the above result, numerical data are pre- 
sented for systems in two spatial dimensions, the first one 
showing an energy threshold, the second one supporting 
breathers of arbitrarily low energy. We consider a Fermi- 
Pasta-Ulam system in two spatial dimensions with inter- 
action potentials of the form -I- and n = 3,4. 
The equations of motion of this Hamiltonian system read 
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FIG. 1: Shape profiles (absolute value of the amplitude) of 
discrete breathers, centered at site 0, in Fermi-Pasta-Ulam 
chains, (a) for a system with phonon spectrum the localiza- 
tion is exponential and decreases when lowering the maximum 
breather amplitude, while (b), for a system without phonon 
spectrum, the localization strength is superexponential and 
remains unchanged. Note the different scales of the ordinates 
(In(-) and ln[ln(-)]). Lines are merely drawn to guide the eye. 



{Xj - X^) [l + \Xj - Xi\' 



(6) 



For n — A, the interaction potential, and hence the 
Hamiltonian 7i, is C°° and we expect the critical dimen- 
sion to be 2, whereas, for rt = 3, we have TL € \ C^, 
and a critical dimension dc = 4 is expected. By numerical 
continuation of periodic orbits from an anti-continuum 
limit [3] J discrete breathers on finite lattices can be com- 
puted numerically up to machine precision. Doing so for 
a set of frequencies approaching the phonon band edge, 
the dependence of the breather energy on its amplitude, 
measured at the site of maximum amplitude, is obtained. 
The result is plotted in Fig. [3 For any finite number 
of lattice sites, a lower bound on the breather energy 
is observed. For n — 3, with increasing system size N 
this bound converges towards zero, whereas for n = 4 it 
converges towards a non-zero value, thus confirming the 
theoretical predictions. 

Systems without a phonon spectrum: A nongeneric, but 
still remarkable, exception from the above reported en- 
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FIG. 2: Breather energy versus amplitude for two-dimensional 
Fermi-Pasta-Ulam systems with periodic boundary condi- 
tions, governed by Eq. JHJ, for the cases n — 3 and n = 4. A 
lower bound on the breather energy is observed. For n — 3, 
with increasing system size A*' this bound converges towards 
zero, whereas for n = 4 it converges towards a non-zero value. 



ergy properties of discrete breathers is observed in sys- 
tems where no quadratic term is present in the inter- 
action potential. Granular materials provide systems 
where such interactions are realized and experimentally 
accessible (see 0] for a review). In such systems, dis- 
crete breathers have been shown to be superexponen- 
tially localized 0, 0, 0, where an upper bound on 
the breather amplitude is given by the inequality < 
aexp(— |i|/6)c('*'/'^^ 0], for some a,b,c,d G M, in one- 
dimensional systems. (The breathers are supposed to be 
centered at site i = 0.) However, the arguments by which 
these results for the spatial decay are obtained should be 
equally applicable in higher spatial dimension, and nu- 
merical computations confirm this reasoning. For these 
sort of systems, the above analysis is not applicable, as 
discrete breathers do not emerge from the bifurcation of 
a band edge plane wave. Due to this, and in contrast to 
the results for systems with phonons reported above, no 
broadening of the discrete breather takes place when low- 
ering its amplitude [Fig.QJb)]. It is this constant localiza- 
tion strength — not the superexponential localization! — 
which allows for the existence of breathers of arbitrar- 
ily low energy regardless of the spatial dimension of the 
system. 
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FIG. 3: Energy versus amplitude for discrete breathers in a 
system without phonon spectrum [Eq. 0] . Independently of 
the spatial dimension d = 1, 2, 3, 4 (from top to bottom), ar- 
bitrarily low breather energies are observed. The data were 
obtained for systems consisting of iV = S'' lattice sites, but, 
due to the strong (superexponential) localization, are indis- 
tinguishable on this scale from those of larger systems. 



This is confirmed numerically for a Hamiltonian sys- 
tem with a harmonic on-site potential and an interaction 



potential of the form 
this system read 



The equations of motion of 



(7) 



In Fig. O for systems of spatial dimension d = 1,2, 3, 4, 
the breather energy is plotted versus its amplitude. In- 
dependently of d, arbitrarily low breather energies are 
observed already in small finite systems. 

Systems without small amplitude breathers: So far we 
have considered systems in which the amplitude of a fam- 
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FIG. 4: Breather amplitude versus frequency for system Q 
with A'^ = 99 lattice sites, free boundary conditions, and 
frequencies close to the phonon band edge (which is at fre- 
quency 2.0). With parameter value \a\ — 0.4 < ^\/3, 
breathers of arbitrarily low amplitude are observed, whereas 
for < \a\ — 1.2 < \/3 the breather amplitude, and hence 
the breather energy, does not go to zero when approaching 
the phonon band edge. 



ily of discrete breathers can be lowered to arbitrarily 
small values, and formerly there has been uttered the 
belief that this were always the case . From a recent 
existence proof for discrete breathers by Aubry, Kopi- 
dakis, and Kadelburg it follows, however, that this is 
not true in general. In a Fermi-Pasta-Ulam chain with 
equations of motion 



3=i±l 



[{xj - Xi) + a{xj - x,)^ + {xj - Xtf] (8) 



and parameter values ja/S < |a| < -s/S, discrete breathers 
do exist, but their maximum amplitudes lie always above 
a certain positive lower bound, even when approaching 
the phonon band edge (Fig. Q . This automatically im- 
plies the existence of a threshold in energy. 
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